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Prompt-photon cross sections in deep inelastic ep scattering were measured with the 
ZEUS detector at HERA using an integrated luminosity of 320pb~ 1 . Measurements 
of differential cross sections are presented for inclusive prompt-photon production as 
a function of Q 2 , x, E^, and rp . Perturbative QCD predictions and Monte Carlo 
predictions are compared to the measurements. 

1 Introduction 

In the study of high energy collisions involving hadrons, events in which isolated high-energy 
photons (often referred to as prompt photons if originating from quarks) are observed provide 
a direct probe of the underlying parton process, since the emission of these photons is un- 
affected by parton hadronisation. Isolated high-energy photon production has been studied 
in hadronic collisions [2 and ep collisions in both the photoproduction [5] and deep inelastic 
scattering (DIS) [1] regimes. The present results are compared to theoretical predictions 
and also to Monte Carlo (MC) models . 

2 Theoretical Predictions 

The cross sections for isolated photon production have been calculated to order 0(a 3 a s ) 
by Gehrmann-De Ridder et al. (GGP) [5]. In the approach of GGP, three contributions 
to the scattering cross-section for ep — ► cyX are considered at leading order (a 3 ) in the 
electromagnetic coupling. One of these contributions comes from the radiation of a photon 
from the quark line (called QQ photons) and a second from the radiation from the lepton line 
(called LL photons). The third contribution is an interference term between photon emission 
from lepton and quark lines, called LQ photons by GGP. In this analysis the contribution of 
the LQ term is reduced to a negligible level by combining e + p and e~p data (the LQ term 
changes sign when e~ is replaced by e + ) and by requiring that the outgoing photon is well 
separated from both outgoing electron and quark. Therefore this analysis took into account 
LL and QQ photons only. 

A calculation based on QED contributions to the parton distributions has been made 
by Martin et al. (MRST) [6]. In the approach of MRST a resummed version of the LL 
contribution is calculated. A photon parton component of the proton arises as a direct 
consequence of including QED corrections in parton distribution functions. This leads to 
a DGLAP-resummed enhancement of the LL prediction relative to that of GGP due to 
the inclusion of QED Compton scattering, 7 p e — > je, where the 7 P is a constituent of the 
proton. A measurement of the isolated high-energy photon production cross section therefore 
provides a direct constraint on the density of the photon within the proton. The non-collinear 
QQ component is not included in the MRST model and so the transverse momentum of the 
scattered electron is expected to balance approximately that of the isolated photon. The 
analysis presented here does not impose this constraint. 
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3 Data samples and Monte Carlo event simulation 

The measurements are based on a data sample corresponding to an integrated luminosity of 
320 ± Spb" 1 , taken between 2003 and 2007 with the ZEUS detector at HERA. The sample 
is a sum of 131 ±3pb _1 of e + p data and 189±5 pb 1 of e~p data with centre-of-mass energy 
y/s = 318 Ge\fl 

DIS events with QQ photon emission were simulated by the MC program Pythia 6.416 
[7]. The LL photons radiated at large angles from the incoming or outgoing electron were 
simulated using the generator Djangoh6 [8], an interface to the MC program Heracles 
4.6.6 [5]; higher-order QCD effects were simulated using the colour dipole model of Ariadne 
4.12[TU]. Hadronisation of the partonic final state was performed by Jetset [TT]. These 
samples were used to study the event-reconstruction efficiency as well as providing the MC 
templates for the signal extraction in Section O and the MC predictions seen in Section [()] 

The NC DIS background used in Section [5] was simulated using Djangoh6, within the 
same framework as the LL events. 

4 Event selection and reconstruction 

The kinematics were reconstructed from the scattered electron^ and the kinematic region 
10 < Q 2 < 350 GeV 2 was chosen. 

Scattered-electron candidates were and required to lie within the polar angle range 
139.8° < e < 171.9°, to ensure that they were well measured in the rear calorimeter 
(RCAL), and have energy, E' e , greater than 10 GeV. To reduce backgrounds from non-ep 
collisions and suppress deeply virtual compton scattering [12j to a negligible level additional 
cuts were made on the reconstructed vertex position, tracks and energy sums. 

Photon candidates were identified as CAL energy flow objects (EFOs [13]) with recon- 
structed transverse energy within the range 4 < E^, < 15 GeV and pseudorapidity within 
the range —0.7 < rf < 0.9. It was also required that at least 90% of the reconstructed en- 
ergy was measured in the barrel electromagnetic calorimeter (BEMC) and that the photon 
candidate was well isolated from reconstructed tracks. 

Jet reconstruction was performed on all EFOs in the event, including the photon candi- 
date, using the kr cluster algorithm [14j in the longitudinally invariant inclusive mode |15j . 
Further isolation was imposed by requiring that the energy of the photon-candidate EFO 
contributed at least 90% of the total energy of the jet to which it was assigned. 

5 Photon signal extraction 

The event sample selected according to the criteria in Section [4] is dominated by background 
events, most notably neutral current (NC) DIS events where a genuine electron candidate is 
found in the RCAL and neutral mesons such as 7r° and r\ decay to photons producing photon- 
candidate EFOs in the BEMC. The photon signal was extracted from the background using 
BEMC energy-cluster shapes. Two shape variables were considered: 

a Hcrcafter 'electron' refers both to electrons and positrons unless specified. 

b Thc definitions of Q 2 , Wx and x as reconstructed using the scattered electron are Q 2 = — (fc — k') 2 , 
W x = (P + k — k' — p-y) 2 and x = Q 2 /(2P ■ (k — k')) where k (fc') is the four-momentum of the incoming 
(outgoing) lepton, p 7 is the four vector of the outgoing photon and P is the four-momentum of the incoming 
proton. 
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• the variable (SZ) = — — w u v e- ' w h erc Zi is the Z position of the centre of the 

ith cell, ^ c i us ter is the centroid of the cluster. Wceii is the width of the cell in the Z 
direction, Ei is the energy recorded in the cell and the sum runs over all BEMC cells 
in the EFO; 

• the ratio of the highest-energy BEMC cell in the EFO to its total EMC energy, /max- 

The distributions of (SZ) and / max in 
data and MC are shown in Fig.Q] The (SZ) 
distribution exhibits a double-peaked struc- 
ture with the first peak at « 0.1, associated 
with the signal, and a second peak at w 0.5, 
dominated by the tt° — > 77 background. 
The / m ax distribution shows a single peak 
at w 0.9 corresponding to the photon sig- 
nal, and presents a "shoulder" extending 
down to « 0.5, which is dominated by the 
hadronic background. 

The number of isolated-photon events 
contributing to Fig. [T] and in each cross- 
section bin was determined by a % 2 fit to the 
(SZ) distribution in the range < (SZ) < 
0.8 using MC-derived templates for the LL 
and QQ signal distributions and the back- 
ground as described in Section [3J In per- 
forming the fit, the LL contribution was 
kept constant at its MC predicted value and 
the other components were varied. The final 
cross sections were derived from the fitted 
events per bin using the MC detailed in Sec- 
tion[3jand the bin-by-bin unfolding method. 

6 Results and Conclusions 

The differential cross sections as a function of Ej., rf , Q 2 and x are shown in Fig. O The 
predictions for the sum of the expected LL contribution from Djangoh6 and a factor of 
approximately 1.6 times the expected QQ contribution from Pythia agree well with the 
measurements, except for an excess at the lowest Q 2 (and associated lowest x). 

The differential cross sections as a function of rf , Ej,, Q 2 and x are compared to the- 
oretical predictions in Fig. [3J The theoretical predictions from GGP describe the shape 
of the E^ and rf distributions well, but the central value of the predictions typically lies 
30% below the measured cross sections. The calculations fail to reproduce the shape in 
Q 2 . As with the MC comparison, the measured cross section is larger than the theoretical 
prediction; this is reflected in an excess of data over theory at low- a;. The MRST predictions 
fail to reproduce the shapes of the measured differential cross sections. However they agree 
with the measurements at large values of Q 2 and x, for backward r\ and for high values of 
Et where the the LL cross section is expected to be a larger fraction of the total. 
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Figure 2: Isolated photon differential cross sections compared to Monte Carlo predictions. 
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Figure 3: Isolated photon differential cross sections compared to theoretical calculations. 
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